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Abstract 

We calculate the amplitude for the forward electroproduction of two light vector mesons in 
next-to-leading order BFKL. This amplitude represents the first next-to-leading order am- 
plitude ever calculated for a collision process between strongly interacting colorless particles. 

1 Introduction 

It is believed that the "gold-plated" measurement for the possible realization of the BFKL 
dynamics [2] is the 7*7* total cross section. However, a prediction for this cross section is not 
yet available with next-to- leading accuracy since the calculation of one necessary ingredient, 
the 7* to 7* impact factor in the next-to-leading order, has not yet been completed after 
year-long efforts [3]. Here we propose as "silver-plated" measurement the differential cross 
section for the 7*7* to two light vector mesons process in the forward case, i.e. for minimum 
momentum transfer. 

In the BFKL approach, both in the leading logarithmic approximation (LLA), which 
means resummation of leading energy logarithms, all terms (asln(s))", and in the next- 
to-leading approximation (NLA), which means resummation of all terms ^^(as In(s))"', the 
(imaginary part of the) amplitude for a large-s hard collision process can be written as the 
convolution of the Green's function of two interacting Reggeized gluons with the impact 
factors of the colliding particles (see, for example. Fig. 1). 

The NLA singlet BFKL Green's function for the forward case of interest here is known 
since several years [4]. Moreover, the impact factor for the transition from a virtual photon 
7* to a light neutral vector meson V = p°, cj, has been recently calculated in the NLA in the 
forward case, up to contributions suppressed as inverse powers of the photon virtuality [5]. 
Therefore we have all is needed to build the NLA amplitude for the 7*7* VV reaction. 

The knowledge of this amplitude is interesting first of all for theoretical reasons, since it 
could shed light on the role and the optimal choice of the energy scales entering the BFKL 

*The content of this contribution is based on Ref. [1], to which we refer for additional details and for a 
more exhaustive list of references. 
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approach. Moreover, it could be used as a test-ground for comparisons with approaches 
different from BFKL, such as DGLAP, and with possible next-to-leading order extensions of 
other approaches, such as color dipole and /cffactorization. But it could be interesting also 
for the possible applications to the phenomenology. Indeed, the calculation of the 7* — > F 
impact factor is the first step towards the application of BFKL approach to the description 
of processes such as the vector meson electroproduction 7*p — > Vp, being carried out at the 
HERA collider, and the production of two mesons in the photon collision, 7*7* — > VV or 
7*7 ^ VJ/"^, which can be studied at high-energy e'^e~ and 67 colliders. 

The process considered here has been studied recently in the Born (2-gluon exchange) 
limit in [6] and in the LLA [7] for arbitrary transverse momentum. In Ref. [7] also an estimate 
of NLA effects has been given in the forward case. 



2 The NLA amplitude 

We consider the production of two hght vector mesons {V = p°, cu, 0) in the collision of two 
virtual photons, 

7*(P) 7*(P') ^ V{p,) V{p2) . (1) 

Here, pi and p2 are taken as Sudakov vectors satisfying pi — p^ — and 2(piP2) = s; the 
virtual photon momenta are instead 

p = api p2 , p = ap2 -pi , (2) 

as a's 

so that the photon virtuahties turn to be = —Ql and (p')^ = —Ql- We consider the 
kinematics when 

s » g2 2 ^ A^^^ ^ (3) 

and 

a = 1 + 91 + o{s-^) , a'^l + ^ + 0{s-') . (4) 
s s 

In this case vector mesons are produced by longitudinally polarized photons in the longitudi- 
nally polarized state [5] . Other helicity amplitudes are power suppressed, with a suppression 
factor ~ mv/Qi,2- We will discuss here the amplitude of the forward scattering, i.e. when 
the transverse momenta of produced V mesons are zero or when the variable t = {pi — p^ 
takes its maximal value to — —Q\Q\/s + 0{s~'^). 

The forward amphtude in the BFKL approach may be presented as follows 

f-^) = / ^*'(^'' / -i^^'^-i- / S (^) « • (^) 

This representation for the amplitude is valid with NLA accuracy. The scale sq is artificial. 
It is introduced in the BFKL approach to perform the Mcllin transform from the s-space 
to the complex angular momentum plane and must disappear in the full expression for the 
amplitude at each fixed order of approximation. Using the result for the meson NLA impact 

factor such cancellation was demonstrated explicitly in Ref. [5] for the process in qiicstion. 

In Eq. (5), $i(gi,So) ^-nd $2(~92i'So) are the impact factors describing the transitions 
7*(p) — >• V{pi) and 7*(p') — >■ V^(p2), respectively. They are presented as an expansion in ag 

^.M = ^1.2 + ')] , ^1,2 = -^¥^v/^c^-i. (6) 
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Figure 1: Schematic representation of the amphtude for the 7*(p)7*(p') — * ^(^1)^(^2) 
forward scattering. 

where ctg = agNc/n, fy is the meson dimensional couphng constant (/p ~ 200 MeV) and Cg 
should be replaced by e/\/2, e/(3\/2) and — e/3 for the case of u and meson production, 
respectively. In the collinear factorization approach the meson transition impact factor is 
given as a convolution of the hard scattering amplitude for the production of a collinear 
quark-antiquark pair with the meson distribution amplitude (DA). The integration variable 
in this convolution is the fraction z of the meson momentum carried by the quark {z = 1 — z 
is the momentum fraction carried by the antiquark): 

Ci'kq') =jdz 011 (^) . (7) 

J q^ + zzQl2 

The NLA correction to the hard scattering amplitude, for a photon with virtuality equal to 
Q^, is defined as follows 



with t{z) given in the Eq. (75) of Ref. [5]. C|^2 (g*^) are given by the previous expression with 
replaced everywhere in the integrand by Qi and Q2, respectively. Below we will use the 

DA in the asymptotic form, = 6z{l — z), both for the simplicity of the presentation and 

because, according to QCD sum rules estimates [8], may be indeed a good approximation 

for the DA of light vector mesons. 

The Green's function in (5) is determined by the BFKL equation 

5^iqi-q2)=ujG^{qi,q2)- J cfqK{qi,q}G^{q,q2), (9) 

where K(qi, ^2) is the BFKL kernel. In the transverse momentum representation (see [1] for 
details), it can be written as 

k = ask'^ + alk^ , (10) 
where is the BFKL kernel in the LLA, k^ represents the NLA correction. 
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The basis of eigenfunctions of the LLA kernel, 



K'\u) = x{y)\y) 



x{v) = 2^{l)-^{- + iv] 



is given by the following set of functions: 

1 



27r2 



iv—iv'—l 



= 5{v - v') . 



(12) 



The action of the full NLA BFKL kernel on these functions may be expressed as follows: 



/5o 



4N^ 



A 



d 



dp 



(13) 



where the first term represents the action of LLA kernel, while the second and the third ones 
stand for the diagonal and the non-diagonal parts of the NLA kernel. We refer to Ref. [1] 
for the expression of x(z^). 

We will need also the |z/) representation for the impact factors, which is defined by the 
following expressions 



du'ci(u'){u'\^ , 



-|-00 

= I dpc2{u){q\p) , (14) 



vrv^ ' -^"^-7 ---- ' ' ^^^^ 

and by similar equations for c{\p) and C2\p) from the NLA corrections to the impact 
factors, dl\q^) and C^^^(g2). 

Using the above formulas one can derive, after some algebra, the following representation 
for the amplitude 



Xnis {A) 
D1D2 



dv 



So 



1 + asiiiR) 



10 dln(^,) 

■ • ^-(^)V21n(/.^) 



du 



We find that 



Tm±tu] 6 



TT 



r[3 ± 2iu] cosh(7rz/) ' 
Using Eq. (16) we construct the series representation for the amplitude 

Q1Q2 TrUsA 1 



D^D2 s 



(16) 



(17) 



(18) 



X 



n=l 



60 + ^oisinRYhn (in ( — ) + dn{sQ,HR)\n-[ — 



So 



So 



n-1 
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where the coefficients 

+00 



/ d^ci(^)c,(^)^, (19) 



Q 



are determined by the kernel and the impact factors in LLA. The expression for the coeffi- 
cients dn can be easily determined from Eq. (16) and is given in Ref. [1]. 

One should stress that both representations of the amplitude (18) and (16) are equivalent 
with NLA accuracy, since they differ only by next-to-NLA (NNLA) terms. 

It can be easily shown that the amplitude (18) is independent in the NLA from the choice 
of energy and strong coupling scales. It is also possible to trace the contributions to each dn 
coefficient coming from the NLA corrections to the BFKL kernel and from the NLA impact 
factors (see Ref. [1] for details). 



3 Numerical results 

In this Section wc present some numerical results for the amplitude given in Eq. (18) for the 
Qi = Q2 = Q kinematics, i.e. in the "pure" BFKL regime. The other interesting regime, 
Qi S> Q2 or vice-versa, where coUinear effects could come heavily into the game, will not 
be considered here. We will emphasize in particular the dependence on the renormalization 
scale hr and sq in the NLA result. 

In all the forthcoming figures the quantity on the vertical axis is the L.H.S. of Eq. (18), 
Ims{A)Q'^/{s DiD2). In the numerical analysis presented below we truncate the series in 
the R.H.S. of Eq. (18) to n = 20, after having verified that this procedure gives a very good 
approximation of the infinite sum for the Y values Y < 10. We use the two-loop running 
couphng corresponding to the value as{Mz) — 0.12. 

We have calculated numerically the 6„ and d„ coefficients for n/ = 5 and Sq — — /ij^, 
getting 



60 = 17.0664 61 = 34.5920 62 = 40.7609 63 = 33.0618 64 = 20.7467 

65 = 10.5698 be = 4.54792 67 = 1.69128 feg = 0.554475 

rfi = -3.71087 (i2 = -11.3057 4 = -23.3879 ^4 = -39.1123 

d5 = -59.207 de = -83.0365 dr = -111.151 dg = -143.06 . 



(20) 



In this case contributions to the dn coefficients originating from the NLA corrections to the 
impact factors are 

rf/^P = -3.71087 4"^ = -8.4361 rff'' = -13.1984 d'^" = -18.0971 

P = -23.0235 d^^ = -27.9877 d^"^ = -32.9676 d^"^ = -37.9618 . ^ ' 

Thus, comparing (20) and (21), we see that the contribution from the kernel starts to be 
larger than the impact factor one only for n > 4. 

These numbers make visible the effect of the NLA corrections: the dn coefficients are 
negative and increasingly large in absolute values as the perturbative order increases. The 
NLA corrections turn to be very large. In this situation the optimization of perturbative 
expansion, in our case the choice of the renormalization scale /ir and of the energy scale 
So, becomes an important issue. Below we will adopt the principle of minimal sensitivity 
(PMS) [9]. Usually PMS is used to fix the value of the renormalization scale for the strong 
coupling. We suggest to use this principle in a broader sense, requiring in our case the 
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Figure 2: Xms{A)Q'^ / {s D1D2) as a function of Iq at fiR = lOQ. The different curves are for 
Y values of 10, 8, 6, 4 and 3. The photon virtuahty Q'^ has been fixed to 24 GeV^ {uf — 5). 

minimal sensitivity of the predictions to the change of both the renormalization and the 
energy scales, /ir and sq. 

We replace for convenience in (18) ln(s/so) with Y — Yq, where Y — ln(s/Q^) and 
Yq — \n(so/Q'^), and study the dependence of the amplitude on 1^0- 

The next two figures illustrate the dependence on these parameters for (^^=24 GeV^ and 
Uf = 5. In Fig. 2 we show the dependence of amplitude on Yq for /ir — lOQ, when Y takes 
the values 10, 8, 6, 4, 3. 

We see that for each Y the amplitude has an extremum in Yq near which it is not sensitive 
to the variation of Yq, or sq. Our choice of /ir for this figure is motivated by the study of 
fiR dependence. In Fig. 3 we present the /ir dependence for Y — 6; the curves from above 
to below are for lo=3, 2, 1, 0. 

Varying /ir and Yq we found for each Y quite large regions in /ir and Yq where the 
amphtude is practically independent on /ir and Yq. We use this value as the NLA result for 
the amplitude at given Y. In Fig. 4 we present the amplitude found in this way as a function 
of Y. The resulting curve is compared with the curve obtained from the LLA prediction 
when the scales are chosen as ^r = lOQ and Yq = 2.2, in order to make the LLA curve 
the closest possible (of course it is not an exact statement) to the NLA one in the given 
interval of Y. The two horizontal lines in Fig. 4 are the Born (2-gluon exchange) predictions 
calculated for /ir — Q and hr — lOQ. 

We stress that one should take with care BFKL predictions for small values of Y, since 
in this region the contributions suppressed by powers of the energy should be taken into 
account. At the lowest order in such contributions are given by diagrams with quark 
exchange in the ^-channel and are proportional in our case to aEMQ^sfv/Q^- -^^ higher 
orders power suppressed contributions contain double logarithms, terms ~ Q;"ln^"s, which 
can lead to a significant enhancement. Such contributions were recently studied for the total 
cross section of 7*7* interactions [10]. 

If the NLA (and LLA) curves in Fig. 4 are compared with the Born (2-gluon exchange) 
results, one can conclude that the summation of BFKL series gives negative contribution to 
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Figure 3: Xms(A)Q'^ /(s D1D2) as a function of Hr at Y=6. The different curves are, from 
above to below, for Yq values of 3, 2, 1 and 0. The photon virtuality has been fixed to 



the Born result for F < 6 if one chooses for the scale of the strong coupling in the Born 
amplitude the value given by the kinematics, /ir = Q. We believe that our calculations 
show that one should at least accept with some caution the results obtained in the Born 
approximation, since they do not give necessarily an estimate of the observable from below. 

Another important lesson from our calculation is the very large scale for (and therefore 
the small itself) we obtain using PMS. It appears to be much bigger than the kinematical 
scale and looks unnatural since there is no other scale for transverse momenta in the problem 
at question except Q. Moreover one can guess that at higher orders the typical transverse 
momenta are even smaller than Q since they "are shared" in the many- loop integrals and 
the strong coupling grows in the infrared. In our opinion the large values of fiR we found 
is not an indication of the appearance of a new scale, but is rather a manifestation of the 
nature of the BFKL series. The fact is that NLA corrections are large and then, necessarily, 
since the exact amplitude should be renorm- and energy scale invariant, the NNLA terms 
should be large and of the opposite sign with respect to the NLA. We guess that if the NNLA 
corrections were known and we would apply PMS to the amplitude constructed as LLA + 
NLA-corrections + NNLA-corrections, we would obtain in such calculation more natural 
values of /ir. 

In the last years strong efforts have been devoted to the improvement of the NLA BFKL 
kernel as a consequence of the analysis of coUinear singularities of the NLA corrections and 
by the account of further collinear terms beyond NLA [11]. This strategy has something in 
common with ours, in the sense that it is also inspired by renormalization-group invariance 
and it also leads to the addition of terms beyond the NLA. These extra-terms are large and 
of opposite sign with respect to the NLA contribution, so that they partially compensate the 
NLA corrections. The findings of the present work suggest, however, that the corrections 
to the impact factors heavily contribute to the NLA amplitude, being even dominating in 
some interval of non-asymptotically high energies. Moreover, by inspection of the structure 
of the amplitude in the regime of strongly asymmetric photon virtualities, one can deduce 
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Figure 4: Ims{A)Q'^ / (s D1D2) as a function of Y for optimal choice of the energy parameters 
Yo and fi^ (curve labeled by "NLA"). The other curves represent the LLA result for Yq = 2.2 
and hr = lOQ and the Born (2-gluon exchange) hmit for Hr = Q and iir = lOQ. The 
photon virtuality Q'^ has been fixed to 24 GeV^ (n/ = 5). 

that also the impact factors generate coUinear terms which add up to those arising from 
the kernel, see e.g. Eqs. (84) and (85) of Ref . [5] . This leads us to the conclusion that in 
the approaches based on kernel improvement the additional information coming from impact 
factors should somehow be taken into account when available. These issues certainly deserve 
further investigation and we believe that useful hints in this direction can be gained from 
the study of the 7*7* — > VV amplitude in the regime of strongly ordered photon virtualities. 

We conclude this Section with a comment on the possible implications of our results for 
mesons electroproduction to the phenomenologically more important case of the 7*7* total 
cross section. By numerical inspection we have found that the ratios 6„/6o we got for the 
meson case agree for n = 1 -j- 10 at 1 -^ 2% accuracy level with the analogous ratios for 
the longitudinal photon case and at 3.5 30% accuracy level with those for the transverse 
photon case. Should this similar behavior persist also in the NLA, our predictions could be 
easily translated to estimates of the 7*7* total cross section. 
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